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We study a self-bound dense object in the hard wall model. We consider a spherically 



o. 

symmetric dense object which is characterized by its radial density distribution and non- uniform 
but spherically symmetric chiral condensate. For this we analytically solve the partial differential 
equations in the hard wall model and read off the radial coordinate dependence of the density 
5_i ■ and chiral condensate according to the AdS/CFT correspondence. We then attempt to describe 

nucleon density profiles of a few nuclei within our framework and observe that the confinement 
scale changes from a free nucleon to a nucleus. We briefly discuss how to include the effect of 
higher dimensional operator into our study. We finally comment on possible extensions of our 
work. 
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1 Introduction 

A self-bound object is a stable bound state in dense matter with strong interactions, i.e., not 
requiring gravity to sustain the bound state unlike a normal neutron star that is gravitationally 
bound. An immediate example of such a state is the nucleus. Nuclei consist of neutrons and 
protons bound by the nuclear force, or quarks bound by the strong interaction. The possibility 
of having strongly bound matter, which would not decay back to ordinary nuclei once formed, 
was first discussed by Bodmer pQ. A well-known example of this kind is strange quark matter, 
in which the energy per baryon number is less than the mass of the nucleon. Strange matter 
has been of great interest ever since the original suggestion that it could be stable [2] and so 
be the ultimate ground state of dense matter. One of the most important consequences of the 
strange matter hypothesis is the possible existence of strange stars, that is compact stars which 
are completely (or almost completely) made of strange quark matter. For details on normal 
neutron stars and self-bound stars, we refer to [3], for example. Another interesting example 
is the kaon-condensed hypernuclei [I], which may be produced in a terrestrial experiment such 
as the heavy-ion collision. The quantities of interest to study the self-bound or gravitationally 
bound object are the structure of the interior and surface of the object, its energy density e(p), 
and pressure P(p), where p denotes the baryon number density. 

In this work, we take a first step towards a realistic description of the self-bound object in 
holographic QCD based on the AdS/CFT correspondence [S]. Recent developments of holo- 



graphic QCD in free space [6] and in dense matter [7] are strongly appealing that holographic 
QCD could be a promising and powerful analytic tool in understanding non-perturbative nature 
of quantum chromodynamics (QCD). It is interesting to notice that holographic neutron stars [8] 
and gravity dual of giant resonances of heavy nuclei [9] are recently investigated. As a first try, 
we consider an oversimplified self-bound object: a structureless and spherically symmetric ob- 
ject. The oversimplified object will be characterized by p(r) and cr(r), where p is the baryon 
number density and a is quark-antiquark (chiral) condensate. Here r is for radial coordinate 
of the spherical object. For the sake of simplicity, we study the simple object in the hard wall 
model |10| [TTj . Since the hard wall model has an infrared (IR) cutoff that is introduced by hand 
in the extra dimension, it is naturally plagued by ambiguities in the IR boundary conditions. 
We will work with several IR boundary conditions in Section [3] and in Appendix [A] Then, 
we introduce a higher dimensional operator |12] into the hard wall model to see the interplay 
between chiral condensate and the density in Appendix [Bj Though it should remain a long way 
to go to a claim that we are describing any realistic physical objects, we comment briefly on 
the nucleon number density profile in nuclei within our framework. Interestingly, we observe 
that the confinement scale changes going from a free nucleon to a nucleus. Our hope is that 
this work would be a good starting point towards a realistic description of the self-bound or 
gravitationally bound objects in nature. 

2 Hard wall model with a higher dimensional operator 

We introduce the hard wall model with U{Nf)i x U{Nf)n gauge symmetry developed in [TUj [TT] 
and its extended version with higher dimensional operators [12]. The action of the hard wall 
model [El EI] is 



S 5 = J d 4 x I ' dzyfg Tr 



--L(Fl + Fl) + \DX\ 2 + ^\X\ 2 



(1) 



where D^X = d^X - iA L ^X + %XA Ril and A L)R = A a LR t a with Tr(t a t b ) = \5 ab . The bulk 
scalar field is defined by X = Xq e 2t7T * , where Xq = (X). The background metric of the model 
is a slice of AdS metric, 

1 



ds = — {r\iLv dx^dx 11 — dz ) , e < z < z m . (2) 



Here 55 is the five dimensional gauge coupling, g\ = I2ir 2 /N c , and e (z m ) is the UV-cutoff 
(IR-cutoff). In [101 [TT] the IR-cutoff z m is fixed by the p- meson mass: l/z m ~ 320 MeV. 

In the standard AdS/CFT approach, the baryon (quark) number chemical potential is in- 
troduced as a background for the time component of bulk U(l) gauge field At [3 [13]. In the 
hard wall model |10[lll|. however, due to specific structures of interaction terms, the role of the 
chemical potential is very much limited. Even though we turn on At to introduce the baryon 
chemical potential, it does not couple to any other bulk fields in the hard wall model. A simple 



way to couple Xq to At is to introduce a higher dimensional operator. As discussed in [12], we 
can add several higher order terms in the hard wall model. Among them, there is only one term 
that could introduce the U(l) chemical potential into the model and that generates coupled 
terms involving At, Xq, and other bulk fields. The term is proportional to XqF v , where Fy is 
the field strength of the bulk U{\) gauge field. Here the U(l) gauge group is the vector subgroup 
of U(Nf)L x U(Nf)ji. Then the modified hard wall model action reads 

S = J d 5 x^gTr [\DX\ 2 + pX\ 2 - -^F v + Cl \X\ 2 F v ^J . (3) 

The equations of motion are 

V 2 X-j^X- Cl XF v = 0, 
V,F^-Ag 2 Cl V,(\X\ 2 F^)=0. (4) 

Since we are going to study a static spherical self-bound object, we take the metric as 

I 2 

ds 2 = — { dt 2 - dr 2 - r 2 (dd 2 + sin 2 6d<j) 2 ) - dz 2 ) . (5) 



z- 



In this background, we solve the equations of motion in Eq. (JU) to calculate the radial coordinate 
r dependence of At and Xq. From At we can read off p(r), baryon number density distribution. 
We can also study the r dependence of the chiral condensate that is basically the boundary value 
of Xq. Expected behavior is that at the center of the object, r = 0, the baryon number density 
will be large, while the chiral condensate will be very small. At the surface of the object, the 
baryon number density will be (much) smaller than the central value, and the chiral condensate 
will be larger than the center value. 

3 Self-bound objects in the hard wall model 

We first consider the case with ci=0. In this case the baryon number density and chiral con- 
densate are quite independent from each other since Xq and A t do not couple to each other, 
which is not the case in any realistic dense matter; the chiral condensate decreases with increas- 
ing baryon number density. However, we can realize this situation in this hard wall model by 
imposing suitable boundary conditions at the center and on the surface of our spherical object. 
This is because the chiral condensate and baryon number density are nothing but integration 
constants of equations of motion for Xq and At. Note that in the hard wall model |10| [TT] 
nonzero value of the chiral condensate was from the boundary condition not from a scalar po- 
tential that supports nonzero vacuum expectation value of X. This feature of the hard wall 
model was improved in |14] , where a scalar potential V(X) is introduced to obtain a nonzero 
value of chiral condensate by minimizing an action. Since we are interested in the baryon num- 
ber density and chiral condensate, we turn on At and consider the vacuum expectation value of 



X. The equations of motion are then 

<9* A t h rcAj h = , 

z r 

2 3^u 3v 2 2<9 r t; 
<9> 1- -^ + <9> H = , (6) 

z z* r 

where Xq(t, z) = v(r,z)/2. Noticing that the equations are separable, we decompose the bulk 
fields as 

A t (r,z) = A{r)B{z), v(r, z) = C(r)D(z) . (7) 

Then we obtain four decoupled second order differential equations that can be solved analytically. 
We first consider the gauge field. The equation for A(r) is given by 

A(r) + -A(r)-l 2 A(r)=0, (8) 

r 

where we use " • " for a derivative with respect to r. Then the solutions read 

A(r) = oo - — , 1 = 0, 
r 

e —lr e lr 

Air) = a 3 h 04 — , I = const. (9) 

r r 



For B{z) we have 



B"(z)--B'(z) + 1 2 B(z) = 0. (10) 



Here " ' " denotes a derivative with respect to z and 

B(z) = b + b -±z 2 , 1 = 0, 

B(z)=b 3 zJ 1 (lz) + b i zY 1 (lz), 1 = const. (11) 

Collecting the solutions, we write the general solution 

*('.^-?)(fc + ^)+£M<fc)(«^+A^) ■ (12) 

The equations of motion for the scalar field and corresponding solutions are 

C(r) + -C(r) - k 2 C(r) = 0, 
r 

C(r) = co--, k = 0, 

r 

e -kr e kr 

C{r) = C3 h C4 — , k = const. , (13) 



and 



3_, ,, / 3 



D{z)" - ^D(z)' + ( -^ + k 2 ) D(z) = , 



z \ z* 

3 



D(z) =d z + diz 6 , k = 0, 

C(z) = d 3 z 2 J 1 (kz) + diZ 2 Yi{kz), k = const. (14) 



The general solution, which is regular at r = 0, takes a form 

v(r,z) = d z + d x z 3 + J2 S * ( Smh r (fcr) ) z 2 J 1 (kz) . (15) 

Expanding the solution near the boundary z = 0, we confirm that the solution follows what it 
should be from the AdS/CFT correspondence as 

v(r, z) ~ d z +L+Y,Sk - S ' mh r {hr) )z 3 + 0(z 5 ). (16) 

Since v or {X) is dual to a quark bi-linear operator of the boundary gauge theory, do corresponds 
to the quark mass up to a constant. For simplicity, we work in the chiral limit, where the quark 
mass is zero. Now we are about to impose boundary conditions to fix integration constants 
above. A well-known plague in the hard wall model is an ambiguity in the IR boundary at 
z = z m . At this moment we do not know how to avoid this problem, and so we make some 
choices of the IR boundary condition below and in Appendix IA1 

3.1 A t (r,zm)=0 

The general solutions in the chiral limit are 

A t (r,z) = (ao-j) (h + |/) + £ zJ x (lz) L^ + A ^\ , (17) 

v(r,z) = cozS + YsSkf^^z'Mkz). (18) 

k ^ r ' 

To ensure the regularity of At at r = 0, we impose a± = and a\ = — /?;. Then the solution 

becomes 

Mr, .)-(/.+ |- 2 ) + E «*«*) (A S -^) ■ ("> 

where [i = aobo, which is identified with the chemical potential of the U(l) charge, and b\ = aob\. 
If we choose an IR boundary condition as At(r, z m ) = and require that two terms on the right 
hand side of Eq. (|19p vanish separately, we obtain 

sinh(j^r) 



A t (r,z) = L- £zA + ^zJ! f^-z) U n 



(20) 



where w n is n-th zero of J\{lz). Now we can read off the baryon number density from this 
solution. At the boundary z = 0, At should be 

A t (r, z) ~ (i + p(r) z 2 + • • • , (21) 

where 



z m „ zz m f \ Z m J 




Figure 1: Density distribution given in Eq. 



Now we impose two boundary conditions at the center r = and at the surface r = R such that 
p(0) = p c and p(R) = pr. Then we have 

2 



Pc 



PR 



P 



+£ 



Pn (Wn 



_P_,ST> / Wn 
?2 + Z> I Q~_ 



/8n 



sinhf^i?) 
R 



(23) 

(24) 



Here /9 C , pr and i? are external parameters whose values are supposed to vary from objects to 
objects. For instance, suppose we are interested in the nucleon density as a function of r, a 
radius of a nucleus which contains A nucleons distributed inside a sphere with a radius R, then 
we can choose 



p c ~0.17fm" 3 , p R ~0 with R ~ 1.2 A 1/3 fm. 
If we consider the lowest mode onhn, n = 1 from (1231) and (1241) , we obtain 



/< 



2 ZmWiPi fwiR 

~PRZ m H — — smh 



2R 

10? I0l 



ft. --=2(p c -p R ) ^-- -^sinhf- 



m 



U>1 



(25) 

(26) 

(27) 



Note that the values of p and /3i are determined by w\, p c , R, pr, and z m . As a result, p is a 
r-independent constant and p(r) is given by 



p(r) = y + 7, srn h ( — r 

Zm 



>2 
-in 



(28) 



In Fig. [H we plot the density as a function of r, where we take pr = 0. 



We do not consider higher modes here since they are expected to have higher energy. 



3.2 v f (r, z r 



Similarly we obtain 



„ (r ,, )=Ei „pii^u M _ z 



(29) 



and 



a(r 



smn | — r 

&Z<mT 



(30) 



Now we impose the following boundary conditions: a{r) = at r = and <r(r) = a® at r = R, 
where a® is the value of the chiral condensate in free space. The assumption for this choice is 
that at the center of a spherical object the baryon number density is high enough to induce 
chiral symmetry restoration and at the surface the density is almost zero. Then we have 



n '"• 

a ° = 2F~H Smh [~ R 



As a simplest nontrivial case we consider the lowest two modes in the Bessel function, 

„2 



x\8i = -x\b 2 , 5 2 = — h5t, 



and end up with 



v(r,z) = 5i 



sinh ( f 1 - r 



xi 



*Jil^z)-^ 



.1 ^ 



J2 I sinh I -^2- r 

•^1 f \ 2m 



z 2 J X ( ^z 



(31) 
(32) 

(33) 



(34) 



As we expand above equation near z = 0, 



air 






Ji3 .2 /sinhf^r 



2l2Zn 



(35) 



where Si is determined from Eq. (|32|) as follows: 



5i = 2ct 



aq 



sinh ^i? 



•H 



sinh ^iT 



i? 



^2^r; 



R 



Now we determined every parameters and a(r) in (j35[) is shown in Fig. [2j 



(36) 



a{r) 
cr 



/? 



Figure 2: Chiral condensate as a function of r 

3.3 Density distribution in nuclei 

Nucleon number density p(r) in a nucleus can be compiled from the charge distribution p c h( r ) 
in nuclei which is measured by elastic electron scattering experiments. With the assumption 
that neutron and proton distributions in nuclei are the same, we can obtain the number density 
profile of a nucleus having A nucleons and Z protons as 



P( r ) = gPch , 



(37) 



where p c h is the charge distribution of a nucleus. A very schematic view of the nucleon density 
profile of several nuclei are shown in Fig. [3j To draw this we use the following formula which is 
valid for A > 20, 



p(r) 



A) 



(38) 



1 _|_ e ( r - R l/2)/a ' 

where po is the normal nuclear density, a = 0.54 fm, and R1/2 = R s — O^A -1 ' 3 fm with 
R s = 1.128A 1 ' 3 fm |15j . Here the radius of a nucleus is given by R ~ r^A 1 ' 3 , where r$ ~ 1.2 fm. 
This relation is easy to understand. Suppose we have A nucleons of radius r inside a nucleus of 
radius R with assumption that both of them are hard spheres, then we can write 



A 



4 s 



(39) 



From this one can have R ~ r$A ' 3 , where tq is a constant. If we assume that nucleons are 
uniformly distributed inside a nucleus, the uniform number density proportional to A/R? is 
independent of the nucleon number A. 

One remarkable point in the figure is that the plateau regime, in which the number density 
is about the normal nuclear matter density, becomes wider with increasing A. 

To show some applicability of our study to a physical system, we explore the nucleon number 
density distribution in nuclei using the result in section 13.11 Since our approach itself is simple 



and we only consider strong interactions, we are not able to explain fine structure of the distri- 
bution, but only a global feature of it. In our approach we have no parameter corresponding 
to the nucleon number A, which may be done with Chern-Simons term. To incorporate the 
A-dependence, we use the relation R ~ r^A 1 ' 3 , and so the size of a nucleus is an input in our 
study. We first show our result with z m fixed in the hard wall model, l/z m ~ 320 MeV, in 
Fig. HI which looks quite different from Fig. [3l We then plot a few density profiles with different 
IR cutoffs in Fig. Here we use a different value of z m from the hard wall model to simulate 
the fact that the plateau in the density profile expands with A in Fig. [5j l/z m ~ 72.8 MeV 
for A = 20, and l/z m ~ 79.0 MeV for A = 50, and l/z m ~ 78.5 MeV for A = 70. In Table 
1, we list several values of A with corresponding l/z m . Here z m is fixed such that for each A, 
our result reproduces the plateau regime of the nucleon number density profile shown in Fig. [3j 
To have a criterion for the closeness of our results to Fig. [31 we consider the root mean square 
(rms) \/ ^~ Jq° \pi( r ) — P2(r)\ 2 dr, where pi(r) is the profile in Eq. (j38j) and P2(r) is our result. 
Here ro is chosen as the largest value of r satisfying pi(r) = P2( r )- For each A we evaluate 
the rms as a function of z m to fix the value of z m that renders the lowest value of the rms. 
We note here that typical errors involved in holographic QCD studies are (10 — 30)%, and so 
we should not seriously take the z m -dependence of A in the table also in Fig. [5J at best, we 
may say that the value of z m for a nucleus can be different that for a free nucleon. As one 
can see from Fig. our results with z m different from the hard wall model is much closer to 
Fig. [3] when we focus on the plateau. If we take our result seriously we can say that varying the 
nucleon number A is related to a change in our IR cutoff z m . Note that l/z m in the hard wall 
model was introduced to ensure the confinement, and so it is intimately related to a confining 
scale. Therefore, we arrive at a qualitative conclusion that the confinement scaleQ ~ l/^m in 
our approach, changes from a nucleon to a nucleus. Though it may not be easy to confirm if 
our finding is true through any experiments, it is interesting to note the claim made to explain 
the European Muon Collaboration (EMC) effect; the quark confinement scale changes in nuclei 
compared to that of free nucleon [16] . 



Table 1 : The nucleon number A and the IR cutoff z r , 



A 


J-/ ■Zm 


20 


72.8 MeV 


30 


77.5 MeV 


50 


79.0 MeV 


70 


78.5 MeV 


100 


77.0 MeV 



2 Here we assume that the confinement scale of mesons and baryons are the same. 




Figure 3: Nucleon density as a function of the distance to the center of the nucleus 

p(r) 



0.15 



0.10 



0.05 




:70 



Figure 4: Nucleon density distribution in nuclei obtained from holographic QCD with z m fixed 
in the hard wall model, where l/z m ~ 320 MeV. 

4 Discussion 

We studied a self-bound dense object in the hard wall model. For simplicity, we considered a 
spherically symmetric dense object which is characterized by its radial density distribution and 
non-uniform chiral condensate. We analytically solve the partial differential equations of motion 
derived in the hard wall model and read off the radial coordinate dependence of the density 
and chiral condensate using the AdS/CFT dictionary. Since we have an ambiguity in the IR 
boundary condition, we tried several different boundary conditions. Though our model study 
might be still far from any realistic physical systems, we tried to understand nucleon density 
profile of a nucleus within our framework. Interestingly, we observed that the confinement scale, 
l/z m in our model, might change from a free nucleon ~ 300 MeV to a nucleus ~ 100 MeV. We 
incorporate the higher dimensional operator into our analysis in Appendix [Bl 

Now we discuss how come the self-bound object forms in the hard wall model. Generally 
speaking, the binding is due to the strong interaction of gluons or gauge fields of a boundary 
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p(l-) 



A = 70 




Figure 5: Nucleoli density as a function of the distance to the center of the nucleus obtained 
from holographic QCD: l/z m ~ 72.8 MeV for A = 20, and l/z m ~ 79.0 MeV for A = 50, and 
l/zm ~ 78.5 MeV for A = 70. 

gauge theory which is encoded in the dual gravity background, the AdS background in our case. 
In our analysis, we did not analyze the stability of the bound object. Therefore it is still not clear 
how we can have such an object in the hard wall model, even though we believe that the binding 
would be supported by the strong interactions or dual gravity background. In fact in our study 
the binding is realized by the boundary conditions in the AdS background. To visualize this, we 
take the chiral condensate in the hard wall model as an example. As pointed out in Ref. [H], no 
scalar potential V(Xq) to support the chiral condensate exists in the hard wall model [101 111] . 
Then, the nonzero chiral condensate was achieved by assigning nonzero value to an integration 
constant which corresponds chiral condensate according to the AdS/CFT dictionary. The value 
of the chiral condensate is determined by observables such as the a± mass or pion decay constant. 
In our case we also make use of the fact that the chemical potential, number density and chiral 
condensate are related to integration constants to assign them suitable values. 

Certainly our present work is leaving many things behind. We didn't properly address the 
issue of stability of our solutions or our bound objects. Therefore there exists a possibility to 
have stable solutions other than the ones obtained in our study. Analysis on the stability of the 
bound object in a well-defined model should follow. The ambiguity in the IR boundary condition 
needs to be removed or to be well under control to obtain any conclusive results. These two 
issues are relegated to our future study in a D3/D7 model [T7]. One immediate generalization 
of our work might be to study nuclear collective excitation modes. Possible connection between 
our bound objects with chiral nuclei in chiral liquids jlS] might be interesting to dig into. One 
may also consider Q-stars [19] or neutron stars; in the case of gravitationally bound neutron 
star one has to introduce the boundary gravity. 
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A More on IR boundary conditions 

In this appendix, we present some more choices of the IR boundary conditions. 
A.l A' t {r,z m ) = 

A' t (r,z m ) = b lZm + Y,Pi (^M.\ {zJi{lz m )) r . (40) 

In order to satisfy a boundary condition A' t (r, z m ) = we set b\ = 0, and {zJ\{y n ))' = where 
Un = lz m - Then we have 

A t (r,t) =/. + £> f ^lM j zJl Q^ . (41) 

Expanding this equation near 2: — » 0, 



/ smh [ -^^-f 
At(r,z)=^+lj2—^\ ^^h 2 + ---, (42) 



2/j r 11 1 



the density can be read off as 

We impose the boundary conditions at r = 0, R as p(0) = p c and p(R) = pr 

Pc=lj2^n, (44) 

sinh ( ^-R ' 






2^Zm r \ R 



If we consider only the lowest mode for n from Eq. (|44[) we have 



.;2 



and pr is determined as 



Pi = 2 Pc ^, (46) 



PR = p^ s[nh (yL R] . (47) 

yiR \z m 
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A.2 v(r,z m ) = 

Now let us consider the scalar solution. At z = z m , we have 

q ^-^ fsmh(kr)\ 2 
v{r, z m ) = c z m + 2^ Ok I J z m Ji(fcz m ) . (48) 

Again we make a simple choice, co = and J\(h n ) = where /i n = kz m , to meet the IR 
boundary condition v(r, z m ) = 

« r ,,) = £ fc p^W£,). (49) 

Near the boundary z — > 0, we can read off the chiral condensate as a function of r, 

/ n lv- r /sinh(f^r)\ h n . . 

2 n \ r J Zm 

B Analysis with higher dimensional operator 

In this Appendix, we take the hard wall model with a higher dimensional operator and present 
basic equations that could be numerically solved. We first turn off the r-dependence to discuss 
the uniform distribution. Then we consider a non-uniform case. Here we will not attempt to 
solve equations of motion numerically partly because due to the ambiguity in the IR boundary 
condition the numerical solution may not mean much. In addition, it is expected that the effect 
of the higher dimensional operator will not change the result obtained in Sec. [3] significantly. 

B.l Homogeneous distribution 

The equations of motion can be written by 

," + i)-^ + ^ + (A-^(^ + i?)), = 0, (51) 



A » i 1 (I i 2c ig lvv> \ (2 2c ig lvv \._ 



(52) 



where " ' " and " ' " denote derivatives with respect to r and z. Since these are not separable, it 
is not easy to solve them analytically. Here we consider only z dependent solution. In this case, 
the equations of motion become simpler as follows: 

A , fl + 2c ig ^W(^ 

\z l-cig£v(zy J 

ii^ 3 ir s ( 3 ^c 1 z 2 A' t (z) 2 \ 
v"{z) - -v'{z) + [^ Z^T^J = ° • ( 54 ) 

From the above, one can easily see that the equation for the gauge field is integrable. Thus we 

can write down the gauge field with v(z) as 

f z x 

A t (z) = fi-p j-TTjrfi , (55) 

Jo l-cig£v(xy 

13 



where the chemical potential \i is given by our IR boundary condition, then the chemical potential 



is 



x 



V = P 1 2 ( \2 dx ■ ( 56 ) 

Jo l-cig£v(xy 



Putting it into Eq. (loll) , we get 

A*) " ~^) + f 4 " tw^tw) W W = ° • ( 57 ) 

z \2r 1/(1 — ci g| ^(z)^/ 

This is the only equation that we have to solve in this case. 

B.2 Perturbation with higher dimensional operator 

Now we assume that the correction from the higher dimensional operators is very small since 
c\ is small [12]. Therefore we can take c\ as an expansion parameter in perturbation. The 
equations of motion Eq. (J3J can be rewritten more explicitly as 

dt + z 3 d± + ^)v = ^v(F^f, (58) 

Of + zd z \d^ A t = -^V, (JF") , (59) 

where we have used a metric ds 2 = - — 2 for convenience. Expanding v and A t as 

v ~ v( ' + civ^ 1 ' + 0(c\) and At ~ Aj. + cij4J + 0(c 2 ) , the corrections can be expressed in 
terms of Green's functions as 

v (1) (z,x) = fdz'dx'G v (z,x;z'x')^v^(F^)) 2 {z',x'), (60) 

J 4z z ^ 

A ( i 1) (z,x) = - J dz'dx'G A {z,x-z , x l )-^V^{v^) 2 F^ t )(z\x ! ) . 

Thus we have to find the Green's functions in order to obtain the corrections. It can be done 
easily by considering the boundary conditions carefully. Assuming a Dirichlet boundary condi- 
tion for At and Neumann boundary condition for v at IR cutoff, we obtain the Green's functions 
as follows: 

00 2 / \ / \ exo I — — — \x — ci? I I 

- ( ,,g,y,g') = -E; » Jife.Wg*') \;- , 1 , ( 6 i) 



n=l 



n * \ ^"m / \ ^m 



^ /( i V rf. Z \ Z^YI J \ Z<TY) r 1 \ X X 



\ exp^ z Jx x\ 

G {z,x;z ,x ) = - > . atA Ji — z Ji — z p — =■ , (62) 

^4^iV^ \z m J \z m J \x-x'\ 

where the normalization constants N% and iVjf are given by 

K = ^ \j[^n? + (l - ^) Jl(x n ) 2 } , (63) 

N* = f^K) 2 . (64) 



Putting these Green's functions into Eq. (|60p . we can compute perturbative corrections to the 
leading order solutions obtained in Sec. El 
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